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The diffraction spectrum of coherent waves scattered from fractal supports is calculated exactly. 
The fractals considered are of the class generated iteratively by successive dilations and translations, 
and include generalizations of the Cantor set and Sierpinski carpet as special cases. Also randomized 
versions of these fractals are treated. The general result is that the diffraction intensities obey a strict 
recursion relation, and become self-affine in the limit of large iteration number, with a self-affinity 
exponent related directly to the fractal dimension of the scattering object. Applications include 
neutron scattering, x-rays, optical diffraction, magnetic resonance imaging, electron diffraction, and 
He scattering, which all display the same universal scaling. 



Scattering is one of the most important methods of observation of structural properties of matter [iJ-pj. Fractals, 
on the other hand, have in recent years received enormous attention as models for the structure of matter J7|—fl0[ . 
Thus the relation between the two is of very general interest, as it provides an essential connection between physical 
observables and the highly intriguing fractal geometry of matter. The first to consider this relation was probably 
Berry JO], who calculated some important averages for scattering by a random self-afhne screen, coining the term 
"diffractal" , for waves that have encountered fractals. In subsequent works, it was demonstrated that diffractals have 
properties that differ significantly from "ordinary" scattered waves. The central feature that distinguishes diffractals 
from ordinary wave fields (where geometrical optics applies), is that they are scattered by objects that have structure 
on all scales, in particular near their wavelength A. This fact gives rise to various scaling laws of the diffraction 
spectrum, reflecting the fractal structure of the scattering object. In contrast, in scattering from crystalline material, 
characteristic interference (Bragg) peaks are observed, which are related through their positions to the underlying 
crystal structure. 

In a large variety of fields scattering processes can be described by a Fourier transform (FT), which relates the 
scattering amplitude to some local density or potential. Examples (to be dealt with in some detail in SecJnJ) include 



scattering |3|, and Helium scattering [[l3|. In all these cases the FT is an approximation, but its generality and 
simplicity have rendered it by far the most widely used approach to scattering problems. The FT is applicable if 
multiple scattering and resonances can be neglected, which is typically the case under conditions of high incidence 
energy. A large bulk of literature, theoretical as well as experimental, exists on scattering in the FT approximation from 
random fractals. The well-known result is that the intensity I(q) decays as a power-law of the momentum transfer q, 
with the exponent related to the fractal dimension D of the scatterer p,|l4]-[l8| . Since they are naturally less abundant, 
much less attention has been devoted to the scattering from fractals which can be constructed by a deterministic set 
of iterative rules. Scattering from such fractals, as well as randomized versions of them, will be the subject of this 
paper. The few examples include Berry's flTl| above-mentioned work; further, mainly in optics, calculations on wave 
transmission Jl^,|2(| and Fraunhofer diffraction f2l|-|24f| , on Cantor-bars, Koch fractals or Sierpinski-carpet like media; 
in x-rays, numerical calculations on scattering by a Menger sponge p5[ , and measurements on diffraction from Cantor 
lattices |2(|. The most extensive treatment is probably due to Allain and Cloitre [p7|-p9|. In Ref. pj|, these authors 
reported an optical diffraction experiment on deterministically generated Cantor bars and Vicsek fractals, and showed 
the resulting structure factor to be self-similar. In Refs. |27],|29|], they analytically solved and discussed properties such 
as band structure and scaling, for the diffraction problem in the case of scattering from a class of fractals similar 
to those to be discussed here. However, their discussion is essentially limited to Fraunhofer diffraction and does not 
include random fractals. Diffractal scattering for probes such as He-scattering or MRI appears not to have been 
discussed in the literature. Thus, there seems to have been no general treatment of the diffractal-FT problem. 

The purpose of the present contribution is to demonstrate that an exact solution for this problem is possible, 
in the case of scattering by objects on an iteratively generated fractal support (see FigB for an illustration of the 
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concept). The class of fractal objects that will be considered here are those that can be generated by a combination 
of dilations and translations. Well-known examples of such objects include (generalizations of) the Cantor set, Vicsek 



fractal, and Sierpinski carpet |3Q]. An operator formalism will be introduced for t his purpose in Sec. [II, which will 



allow the treatment of diverse scattering conditions. It will appear in Sees. IV -VII that whereas some details of the 



diffraction spectrum are context-sensitive (i.e., determined by a form- factor), the ove rall structure is determined by a 



universal, context-independent scaling relation. This conclusion is un altered (Sec.|VTj) by the introduction of a fractal 



dimension-preserving randomness. Following this finding, Sees. VIII IX attempt to address the connection between 
the results derived here for a somewhat artificial class of fractals, and the standard (power-law) expressions used 
to interpret scattering data from natural fractals, such as self-similarity displaying porous materials, aggregates, or 
ramified structures. 

The fractal scattering object can be generated in two essentially different ways: from bottom-up (henceforth BU 
- iterative inflation) or from top-down (TD - iterative refinement). In the former, the smallest unit remains fixed, 
and structure appears at ever larger scales, limited of course by a natural upper cut-off. This structure is reflected 
at ever smaller scales in momentum space. A fixed point is reached where further spectral details are indiscernible, 
either due to experimental resolution constraints or when the wavelength becomes larger than the upper cut-off. 
To every iteration of the fractal support there corresponds a diffraction spectrum. Subsequent diffraction-spectrum 
iterations may be equated when the fixed point is reached. The BU description is appropriate, e.g., in the case 
of a fractal formed around a single nucleation center in a deposition process, as coverage is increased. In the TD 
case, the total system size is fixed and structure appears at progressively smaller scales, limited by a natural lower 
cut-off. This description is probably more appropriate for the physical formation of fractal structures by removal 
of material (pore-fractals (31]). If the wavelength A of the incident waves is fixed, there will necessarily be another 
fixed point, where structure develops below A, and further fractal details are indiscernible. Another possibility is that 
A becomes smaller than the lower cut-off. In both cases subsequent iterations of the diffraction-spectrum can then 
again be equated. Consequently, in both BU and TD cases, as will be shown here, the diffraction spectrum becomes 
(approximately) self-affine, and the self-affinity (or Holder) exponent is simply related to the fractal dimension of 
the scattering object. The central new result derived here is that this conclusion is unaltered neither by the physical 
identity of many scattering probes, nor by the introduction of a fractal dimension-preserving randomness. 

II. FOURIER-TRANSFORM RELATIONS FOR COHERENT WAVE SCATTERING 

The purpose of this section is to summarize the relation between the structural properties of the scattering set and 
the observable diffraction spectrum, for various physical examples to which the FT is applicable. The ultimate goal 
is to show that in spite of the apparently very different way in which the interaction (potential) between the wave 
and the scattering object enters the formulation in each of the cases considered, there are certain universal features in 
the scattering intensities, which reflect only the underlying fractal geometry of the scatterer. The different examples 
are presented b elow in increasing order of computational probe-object interaction complexity. Thus, whereas ne utron 
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scattering (Sec. [I A) involves merely a discrete Fourier sum over the nuclear coordinates, electron scattering (Sec 
requires the FT of a potential which is a functional of the local electron density, and He scattering (Sec. 
necessitates the Fourier transformation of a functional of the interaction potential itself. Yet, it should be emphasized 
that the results presented in this work, all pertain exclusively to local interaction potentials (as holds for all the 
examples considered below). Non-local potentials have been successfully considered in the literature as well, mainly 
in low-energy nuclear problems [ ^2|j33| . 

Notation and conventions: The momentum transfer is denoted by q; the scattering amplitude by f{q)\ spatial 
vectors by f = (x,y,z). Elastic scattering is assumed throughout. As emphasized in each of the subsections below, 
the FT is essentially always the consequence of a high-energy approximation. 

A. Neutron Scattering 

Neutrons may couple by virtue of their spin to magnetic moments. However, the interaction of interest in the 
present context, i.e, which gives rise to a Fourier integral, is with non-magnetic material, where neutrons are scattered 
by the nuclei. Due to the extremely short range of the strong force, t his p rocess is treated almost exactly in the Born 
approximation. The neutron-nucleus interaction potential [see Eq. ( |2.7| )] is essentially a delta function (the "Fermi 
pseudo-potential" g), so that if the nuclear positions are {n}, then 
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The accuracy of this expression depends on the extent to which one may neglect incoherent scattering due to isotopes, 
and inelastic diffraction due to variation of the structure with time (thermal vibrations or atom diffusion). 



B. X-Rays 



The well known Laue derivation |34|], yields the relation 

f(q) = J drn(r)e^ p (2.2) 

between the local electron concentration n[f) and the x-ray scattering amplitude. The assumptions underlying 
the Laue derivation are essentially that the polarization and electric field intensity are linearly and locally related 
by the dielectric susceptibility which itself is frequency-independent. Furthermore, at the inherently high x- 

ray frequencies, x <C 1, which allows for a decoupling of the equations resulting from the attempt to solve the 



electromagnetic wave equation in the crystal lattice, and yields Eq. (2.2) 



C. Optical Diffraction 



The FT arises in optics in the case of Fraunhofer diffraction. This holds when both source and observation point are 
located very far from the aperture, although some more general conditions exist Q). The Fraunhofer formula results 



from the small-wavelength Kirchhoff theory 35 1, which solves the wave-equation under Huygens-Fresnel boundary 
conditions. The assumed smallness of the optical wavelength in comparison with the dimensions of the diffracting 
obstacles implies that in optical diffraction, the BU fractal construction is more natural. Essentially, Fraunhofer 
diffraction occurs when a coherent light wave is scattered by an object with transmission function t(r), and the light 
amplitude is obtained by a coherent superposition 

f(q) = C I drt(r)e 1 ^. (2.3) 



D. Magnetic Resonance Imaging (MRI) 



Suppose the local nuclear spin density in a sample is p(r), and that an oscillating magnetic field with local Larmor 
frequency w(r) is applied to it. It is conventionally assumed in MRI that the Larmor frequency is linear in the nuclear 
spin coordinates: 

cj(rO= 7 |-Bo|+7G-r, (2.4) 

where 7 is the gyromagnetic ratio, and B is the polarizing field, much larger than the linearly varying gradient field, 
of which G is the gradient. In practice, heterodyne mixing eliminates the term 7I-B0I) and the integrated MRI signal 
amplitude can be written as 

f(t) = J drp(r)e l -^ pt . (2.5) 

A reciprocal space vector q = "fGt is introduced |l2] ] , showing that g-space may be traversed by moving either in time 
or in gradient magnitude, so that 



(2.6) 



Eq. (2.6) assumes rapid signal acquisition (after the excitation pulse), so that spin relaxation, dipolar and scalar 
coupling, or spin translation, do not distort the signal. 
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E. Scattering of Electrons from Atoms 



Here one often applies the Born approximation, 



in 



f(q) = -^J dre~^ r V(r), (2.7) 

valid at high energies and assuming that the electron (of mass m) sees a fixed electrostatic potential due to a charge 
density n(r), 

V(r) = -ej & ^L. (2.8) 
This expression neglects the possible polarization of the atom by the incident electron, as well as exchange effects ||. 

F. He Scattering 

The He-surface scattering problem has been successfully treated within the Sudden approximation |j(|[j7|], which 
assumes a high perpendicular momentum change compared to the momentum change parallel to the surface (essen- 
tially a high energy approximation). Under the presence of an arbitrary He-surface potential U{R,z), the Sudden 
approximation yields the scattering amplitudes as 

f(Q) = jJ dRe^e 2 ^, (2.9) 
where the phase-shift function is given in the WKB approximation by 



/>oo 

r)(R) = dz 

h{R) 



2m 



k l - -^U(R,z)) -k z 



1/2 



k t £(R). (2.10) 



Here R = (x, y), Q = (q x , q y ), and k z is the wavenumber component normal to the surface. The turning points £(i?) 
are obtained as solutions to the energy equation 

U[R,^R)} = ^, (2.11) 

with m the mass of the He atom. Effects such as resonances, multiple collisions and dynamic polarization are neglected. 
For a hard-wall potential 



U(R, z) 



: z > £(R) 
co : z < 



so that from Eq. ( 2.10| ) it follows that in this case: 

r)(R) = -k z £(R), (2.12) 

as in the eikonal approximation in optics. 

III. GENERATION OF FUNCTIONS ON FRACTAL SETS BY DILATION AND TRANSLATION 

OPERATORS 

Having seen the generality of the FT in diffraction problems, the generation of the scattering fractal support is 
given next. The construction to be described below is in the spirit of the Iterated Function System formalism of 
Barnsley S. 
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A. Simple Example 



Consider first as an introductory example the construction of a characteristic function on the usual (ternary) Cantor 
set (Fig.|l|, left): One first contracts the generator (zero-order iteration), £o(x) = I (0 < x < L), by a factor 3, and 
then places one copy of the contracted version at the origin, and another translated by 2L/3 from the origin. This 
can be generalized to contractions by a factor 1/s (0 < s < 1) and translations by aL. The corresponding TD fractal 



operator is (the reason for using the adjoint will become clear in Sec. IV): 



where the translation operator is defined as 



and the dilation operator is defined as 



^ = (l + T_ a )C 1/ „ (3.1) 
T a h{x)=h(x + aL), (3.2) 
C s h(x) = h(sx). (3-3) 



T a shifts the function it operates on by an amount aL to the left, and C s stretches the function by a factor of 1/s. 
When applied in the inverse sense as required by the definition of , it is easily checked that £i fx) = £o (x) — 
6o (x/s)+ fo [fa: — aL) / s] , and that £ n (x) = (T^) n £q{x) is indeed an rt th iteration stepped Cantor surface, as illustrated 
in Fig|j]. Barnsley |3^] and Vicsek |59[ provide a general theorem for the calculation of the fractal dimension D of 
such iteratively constructed fractals; D is the solution of the equation 



£■ 



*? = 1 (3-4) 



where Sj are all the contraction factors. Thus in the present case: 
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i=l v ' ' 

To derive the algebraic properties of the above operators, it is convenient to express them in exponential form. T a 
has the well known momentum-operator representation 

T a =e aLd *. (3.6) 

This can be used to find a similar representation for C s : Let fi = ln(s), y — ln(x), and g(y) = h(x). The argument of 
h(s x) can then be expressed in terms of a sum: h(sx) = h[exp(y + /x)] = g(y + /i). But this is exactly in the form of a 
translation, so that using the representation of T a one finds: g(y + fj) = exp(fid v )g(y). Noting that d y — di n r x ) = xd x , 
one obtains the desired representation: 

C, = e M«)««.. (3.7) 
From here, using d x = — d x and d x x = 1 + xd x , it is easily seen that 

# = T_„ 

C\ = -C 1/s . (3.8) 



B. General Construction of Functions on Fractals 



The above formalism for TD fractals can easily be extended to arbitrary dimension, as well as to BU fractals Let 
r = (xi, .., Xd) be a vector in d dimensions. Then the generalization of the ID translation and dilation operators is 

T s h(r) = h(r + aL) 

C s h{r) = h{sr). (3.9) 
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In exponential representation, it is easily seen that: 

T s = e LSV 

Cs = £ >«?-v ( 310 ) 
A very wide class of fractals can be generated by using a single contraction factor s JiCj : 

k 

^ = (1 + £t_3,)C v> . (3.11) 
i=l 

For example, the Vicsek fractal (39|, results by choosing 
a = 1/3; {Si} = {(2/3, 0), (1/3, 1/3), (0, 2/3), (2/3, 2/3)}, 
whereas the Sierpinski carpet is generated by 

s = 1/3; {a,} = {(1/3, 0), (2/3, 0), (0, 1/3), (0, 2/3), (1/3, 2/3), (2/3, 1/3), (2/3, 2/3)} (Fig.§. Eq. @ for the calcu- 
lation of the fractal dimension applies again, and one obtains in the present case: 



The BU fractal is most easily derived by employing the general fractal operator [Eq. ( (3.11 )1, and the observation 



that repeatedly expanding the TD fractal achieves the desired result. Thus the general BU fractal-operator is: 

Gt = {C s ) n {^) n . (3.13) 

Note that with this definition, it is guaranteed that the smallest building-block making up the fractal is of unit length. 
Since the expansion is one-sided, the fractal thus obtained is semi-infinite. 

For future reference it is convenient to note, using Eqs.( |3.8| ) for Tl and C s , that: 

k 



i=l 

\/s 



g n = s- dn r i c[\ s . (3.i4) 



IV. INTRODUCTORY EXAMPLE: ID, HARD- WALL He SCATTERING FROM A CANTOR SET 

With the fractal operators defined, a simple, but prototypical diffractal-FT problem can now be discussed. One 



may, e.g., consider ID He scattering in the presence of a hard-wall potential [Eqs.( pT9| ),(2.12)], with the shape-function 



: else. 



C n denotes the n approximation to the Cantor set. 



A. Calculation of the Intensity Distribution 

1. TD Case 

Denoting the phase-shift of a He atom with perpendicular wavenumber k z and striking a step of height h, by 

$ = -2M, (4.1) 

one notes that exp(iQt; n (x)/l) — exp(z<i>) for x€C n and 1 otherwise. This calls for a normalized chara cteri stic function 
on the Cantor set. Such a function is just j(J^) n £o(%)- Therefore the scattering amplitude is (Eq. ( |2.9| ) in ID), 



G 



Uq) = J,[ dx &iqx ^T^ n ^ + dxe iqx (i ~ yK-^r&O*)]) 

e 4 * - 1 r L , r „„ , . 1 ri 



cfe [re^]6(i) + - dx e«*to(x). (4.2) 



LI J J ^ uv ' L 

The last term is evidently just the specular contribution, and will henceforth be a ssum ed subtracted out. The 
penultimate term contains the fractal operator, which in the present case equals [Eq. ( fig )] 

T = sC s (t + T a ). (4.3) 

What remains is to calculate T n e lqx : 

Te lqx = sC s [(l + e i9Qi ) e lqx ] = s (l + e i<?ai ) e i<?s:E 

J^V 9 * _ s ^ 1 + e igai) :Fe ig S x = fi 2 ^ + gigai^ ( 1 + £ iq s ah ) & iqs 2 x^ ^ 

from which the general pattern can be inferred: 

n 

■pnjqx =s n e iqs"xY[ f 1 + £ iqaLsi-^ _ (45) 
J=l 

This prototypical expression, or slight variants of it, will appear repeatedly when more complicated cases are treated 
in later sections. Before the intensities are obtained, the question of normalization must be addressed. Since the 
Cantor set and its generalizations discussed here have measure zero, the intensity is expected to vanish. This can be 
avoided if the intensity is normalized to the relative length occupied by the Cantor set support at the n th iteration. 
There are 2" steps in the set, each of length s n L, resulting in a normalization factor of L(2s) n /L. 

The integration leading to the scattering amplitude [Eq. (4.2)] can now be performed, yielding, after normalization: 



^ = (2^^ ( e ^ L - n (i + * iqaLsi ) ■ ^ 

The last result bears some resemblance to the (complex-) Weierstrass-Mandelbrot function 0, 

oo 

W(q) = (1 - w 2 )~ 1/2 J2 w ° ( e2 " S3 q ^ 1 )-' ( 4J ) 
j=— oo 

which suggests that the off-specular amplitude, as well as the intensity, 

In(q) = \!n{q)\\ (4.8) 
may be self-affine functions. Before this is investigated, consider first the BU construction. 



2. BU Case 



Essentially, all that needs to be do ne is to replace the TD operator T n in the previous subsection, everywhere by 
the BU operator Q n . From Eq. (3.14) this operator is in the ID case: 



When this is applied to the Fourier basis-set one finds 

g n e iqx = s - n r l e ls ~ nqx = 



n ' 



s J 71 1 q aL 



(4.9) 



(4.10) 



where the last equality follows from the general result for J^ n e z i x [Eq. ( |4.5[ )]. As for normalization, since the fractal 
grows indefinitely in the BU case, it is most convenient to normalize the intensity by the number of elementary units. 
This is 2 n for the n th iteration. 
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In 



anticipation 



of 



the 



more general treatment of Sec. VI, the scattering amplitude f n (q) = j; J dx exp(iqx) exp[— 2i k z £, n (x)] can now 
be written as 



/«(?) = 



I n_1 



2"/ 



F(q) = - j dxe^Mx) 

(j) {x) = e' 2lk ^ oi - x \ (4.11) 
where F(q) can be interpreted as a /orm factor and the term in square brackets as a structure factor S(q) EJ . 



B. Recursion Relation and Self- Affinity of the Off-Specular Intensity Distribution 

Similarly to the fractal sets described above, self-affine functions can be constructed iteratively, for example as 
deterministic models of random walks |3l|j39| , ^2]| . At each stage, a function of this type satisfies the recursive scaling 
relation 

h n+1 (x) = b- a h n (bx), (4.12) 

and becomes rigorously self-affine in the limit n — » oo. a is denoted the Holder, or self-affine exponent [Q. An 
analogous recursion relation will now be derived for the off-specular amplitudes and intensities I n (q) of the previous 
subsections. In the n— >oo limit, these are therefore also self-affine functions. 



1. TD Case 



Using the result derived previously for the scattering amplitude [Eq. (4.6)1, the intensity satisfies 



In{q) 



(2sY 



2 

qL 



[1 -cos($)][l - cos (qs n L)] 



n-l 

2" JJ [l + cos (gaLs*)] 

3=0 



(4.13) 



The recursion-scaling relation follows once it is recognized that the scale factor b from Eq. (4.12) is the dilation factor 
s in the present case: 



l(q) = sl n (sq)[l + cos(qaL)}. 



(4.14) 



Clearly, due to the presence of the cosine factor, this is not in the form of the self-affine scaling relation of Eq. (4.12) 



where a constant factor multiplies the n iteration. However, in the TD case, successive fractal iterations will result 
in successive diffraction spectra that differ at ever larger q scales. q ma x, the largest possible q, is fixed by energy 
conservation, irrespective of the structure of the scattering fractal set. Therefore, when the finest fractal detail, Ax n , 
becomes smaller th an 2it /q maxi it becomes physically reasonable to equate successive iterations. For these to match 
in the sense of Eq. ( 4.12 ), the simplest criterion is to require equality of th e int ensities in the vicinity of the specular, 
q — > (at the price of mismatch increasing with q). Proceeding thus, Eq. (4.14) will be in the form of the self-affinity 
relation [Eq. ( 4.1 2| ) ] if cos(<7aL) is evaluated at q = 0. For then one finds 

I n+1 (q) tts- ai I n (sq) (4.15) 



where 



Oil 



D - 1 



(4.16) 



with D the fractal dimension of the Cantor set, Eq. (3.5). Thus, the self-affinity exponent of the intensity spectrum 
is related to the fractal dimensi on of the object scattered from. The r eason for the specific form of the expression 

(4.15) produces the required scaling can be 



The accuracy with which Eq. 



for ax will become clear in Sec. VI B . 

Fig.|. Plotted there are the intensities for He scattering from a ID, hard- wall step function arrangement on 

Significantly, the intensities of all maxima (not just the 
. This situation can only be expected to improve as n is 



seen m 

two different Cantor set supports (see caption for details) 
specular, corresponding to q — 0) are accurately reproduced 
increased, demonstrating the self-affinity of the spectrum. 
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2. BU Case 



The scaling relation in this c ase is somewhat different from the TD fractal. From the scattering amplitude calculated 
for the BU fractals [Eqs.( |4.1l|) 1 one finds 



— -— = l + cos(s qaL)\ (4.17) 
Ia(q) Io{q/s) 

implying that scaling is obeyed to within the form-factor (i.e, only the structure- factor, not the intensity, is fully 
scale-invariant). For a BU fractal, features in successive diffraction spectra develop at ever smaller q scales. Beyond 
the experimental q-space resolution, it is physically reasonable, as in the TD case, to compare successive iterations, 
and to require the intensities in the vicinity of the specular (q — > 0) to be equal. Substituting 1 for cos(q aL), it is 
now found that 

low io{q/s) 

where 

a 2 = D (4.19) 

with D again the fractal dimension of the Cantor set, Eq. ( |3.5[ ). This is demonstrated in Fig.[|, where the scaling 
recipe with a 2 is seen to hold with high accuracy. 



3. Numerical Check of the Self- Affinity 



To further test the self-affinity, the Holder exponents of the structure factors for He scattering from hard-wall step 
functions on an n = 8 ternary and n — 6, s — 1/5, a — 4/7 Cantor set (Fig^) were calculated, using the reliable 
and accurate epsilon-variation method The result is shown in FigJ^. The respective regression-slopes of 1.367 

and 1.59, from which self-af finity exponents of 2 — 1.367 = .633 and 2 — 1.59 = .41 are obtained, compare favorably 
with the prediction of Eq. ( |4.19| ), yielding a 2 = ln(2)/ln(3) = 0.631 and a 2 = ln(2)/ln(5) = 0.43. Significantly, 
the log-log plots are straight over two orders of magnitude, and the higher order iteration (n = 8) yields a more 
accurate exponent. Note further that this experimental- like analysis yields the same self-affinity exponent as the 



scaling analysis leading to Eq. (4.19), without the q = approximation. 

In conclusion of this section, the analysis of both TD and BU fractals suggests that (1) the scattering intensity from 
a fractal surface is itself (approximately) self-affine, and (2) the fractal dimension of the scattering surface manifests 
itself simply through the Holder exponent of the scattering intensity. Hence an analysis of the scaling properties of the 
scattering intensity should reveal if the scattering surface is fractal, and if so, what its fractal dimension is. 

In the following sections it will be shown that this conclusion holds for the general (arbitrary local potential, any 
dimension) diffractal-FT problem. However, first a commutation property of the operators under discussion must be 
established. This property will make it possible to demonstrate that the scaling discussed above is indeed independent 
of the nature of the scattering probe, and is instead exclusively determined by the geometry of the fractal scattering 
object. 



V. COMMUTATION PROPERTY OF CHANGE-OF- VARIABLE OPERATORS 

The translation and dilation operators T a and C s can clearly be regarded from their definition [Eqs.( |3~2] ),([3~5|)], as 
" change- of-variable operators" (COVO). Let V be a general COVO, i.e., 

Vp = V ' : j/(a;)=J#(aO] (5.1) 
The purpose of this short technical section is to prove that the following commutation relation holds for COVO: 

nfW)\} = f\P9){x)] (5-2) 

To prove this, consider the LHS: Let 







f[g(x)} = ho(x), (5.3) 
and note that the LHS is Vh = h' . But, by Eqs.(U),(f]|, 

h' (x) = ho[4>(x)] = f{g[<t>(x)}} = f[g'(x)\ = f[(Vg){x)\, (5.4) 



which is identical to the RHS of Eq. (5.2), so that the commutation property holds. 

Note also that nothing in the above discussion restricted the result to ID: the commutation property holds in arbi- 
trary dimension. Thus a useful corollary follows immediately. Let dG(x, z)/dz = g(x, z). Then by the commutation 
property: 

V ( <(X) g(x, z)dz = V{G[x, {(x)\ - G[x, £(x)]} = G{</>(x),S[4>(x)]} - G{c(>(x)^[cf>(x)}}; 

r VQ(x) rCMx)] 

Vg(x, z)dz = / g(<f>(x),z)dz = G{<t>(x),([cf>(x)]} - G{<f>(x),Z[cf>(x)]}, 



so that: 



V / g(x,z)dz = / Vg(x,z)dz. (5-5) 
JtW ' Jvt{x) 



VI. GENERAL DETERMINISTIC DIFFRACTAL-FOURIER TRANSFORM PROBLEM 

The tools are now prepared to consider the general, deterministic, diffractal-FT problem. This will require the use 
of the commutation property of change-of-variable operators to treat the variety of scattering probes discussed in 
Sec.[n] in a unified way. 

A. Structure of the Fourier Integral 

The ingredients entering the general problem are (1) the FT relations from Sec.||, and (2) the realization that in 
every such instance, the fractal structure may be introduced into the problem by the repeated application of fractal 
operators to a generator £,a(f). The structure of the general scattering amplitude is therefore: 

f n (q) = J dre^Un(r), (6.1) 

where quite generally (and similarly for the TD case with T replacing Q), 

Mr) = A[GiUr)i (6.2) 



with A an operator to be specified next. For example, in the relatively simple x-ray case [Eq. ( |2.2| )1, A is the identity; 
= n o(r) is the zero-order local electron den sity ; and <^ n (^0 is the electron density on the n th iteration of the 



fractal support. The electron scattering case [Eqs.(2.7),([2.q)] is more complicated, since now A is an integral operator 



acting on Gl l {no(r l )/\r — r"|). But the COVO commutation property, in particular Eq. (5.5) with consta nt i nt egrat ion 
limits, ensures that Q\ can be taken out and put in front of the integral. The He-scattering case [Eqs.(2.£),(2.10)] is 
probably the most complicated, since there the fractal operator G acts at several places simultaneously and A is an 
integral operator with a functional limit. Nevertheless, the COVO commutation property and its corollary simplify 
the problem to the extent that G^ may be pulled out again: 

/n (Q) = JL y d i[ e iRQ e 2ivn{R) = jjdR (G n e^ M ) e 2l ^\ (6.3) 

Indeed, it should now be evident that this is the general structure of the (local-potential) diffractal-FT problem: the 
fractal operator can always be moved from the context-specific part (£o) to operate on the Fourier basis-set, so that 
generically 
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fn(q) = J df (g n e if *) ^ (r). 



(6.4) 



This is the general structure of the Fourier integral: a context-specific part embodied in the integrand of the form- 
factor, 4>o(r), and a generic part common to all diffractal-FT problems, found in the operation of the fractal operator 
on the Fourier basis-set. What remains, in order to understand the universal scaling behavior, is to investigate this 
latter part. 



B. Scaling of the Fourier Integral 



The fractal operators to be dealt with here are of the general form given in the TD and BU cases (Eqs.(3Tl),Eq 



(3.13) respectively). The ensuing analysis closely follows along the lines of the simple, ID case treated in Sec. IV 



1. TD Case 



Repeating the ID calculations of Eq. (4.4), one finds that now the application of the fractal operator yields: 



n-l / k \ 
j=0 \ i=l J 



(6.5) 



The normalization in the general case is to the relative volume occupied by the fractal, i.e., (s n L) d (k + l) n /L d , since 
there are (k + l) n fractal elements at the n th iteration, each with volume (s n L) d . Introducing a form factor, 



F(q) = J dfe l ^M^), 
the scattering amplitude assumes the following generic form: 

n-l / fe 



fn(q) 



id 



(s nd (k + i) n y/ 2 ' 



n l + J2e ls3 ^ S ' L F{qs n ) 



(6.6) 



(6.7) 



3=0 



i=l 



Thus the exact scaling relation for the intensities reads 



ln+i(q) 



k + 1 



1=1 



In(qs), 



(6.8) 



and by employing the recipe used and justified in Sec. IV for ID, of evaluating the exponential terms at q = 0, one 
obtains the approximate self-affinity relation 



In+M W (k + l)s d I n (qs). 



(6.9) 



Expressing this through the Holder exponent as in the ID case, J n -|_i(q') w s _Ql I n (sq), on e find that the universal 
relation between the self-affinity of the intensity spectrum and the fractal dimension [Eq. ( |3.12| )], for a TD fractal 
support, is: 



a± = D — d, 



(6.10) 



in agreement with the ID case. 



11 



2. BU Case 



In strict analogy to the results in ID [Eqs. ( 4. 1C ), ( 4.11 )1, one finds in the d-dimensional BU case: 



g n e iqr = s- dn T n e ls q - r 



n i+E' 

j=i \ i=i 



iq-r 



(6.11) 



(fc + I)™/ 2 



n h+E- 



J'=0 



t is J n q-dtL 



(6.12) 



dr i 



r >o(r). 



(6.13) 



The normalization reflects that there are now (k+ l) n elementary units at the n th iteration. Consequently, the scaling 
of the intensities is: 



In+i{q) 



1 



k + 1 



;.=i 



where again 



In(q/s) _ , , jnCg/f) = -q 2 Jn(g/s) 



a 2 = D 



(6.14) 



(6.15) 



with 13 the fractal dimension [Eq. (|3.12|) 1 of the BU fractal support. Interestingly, it thus appears that the embedding 
space dimension does not enter the scaling in the BU case. This fact remains to be explained on physical grounds. To 
visualize the features of the intensity distribution in this case, Fig. [| displays ID sections of the results of He scattering 
calculations from Ag adatoms centered on a sixth generation Sierpins ki ca rpet with an underlying Pt(lll) s urface 
(BU version of FigJ|). These results were obtained by employing Eq. ( 6.1 2| ) for the structure factor, and Eq. ( 6.13h 
for the form factor. The latter was calculated in the Sudden approximation with a realistic He/Ag/Pt potential |45|] 
for cp (R) = exp[2ir](R)]. 

To summarize, it was shown that irrespectively of the nature of the probe, for the scattering of a coherent wave by 
a deterministic fractal support, the intensity spectrum is approximately self-affine, with a Holder exponent trivially 
related to the fractal dimension of the support. 

The next generalization, necessary to approach realistic situations, concerns the effect of randomness. 



VII. SCATTERING FROM A RANDOMIZED FRACTAL SUPPORT 

Realistic fractals always contain some element of randomness 46[|. For example, in DLA |l7| the adsorbing particles 
perform a random walk and the resulting fractal is consequently random. Thus it is of major interest to introduce 
some randomness into the fractals under consideration, and to investigate its effect on the conclusions reached so far 
regarding the scaling properties of the intensity distribution. In order to meaningfully introduce randomness, it is 
useful to preserve the FD of the support. Otherwise the fractal dimension is not a useful descriptor of the scattering 
object. This preservation of the fractal dimension can be achieved by keeping the constant, single contraction factor, 
but allowing for a distribution of translations. The translations will be chosen independently from a given, but 
arbitrary, probability distribution -P(a), with normalization 

k 

l[d4p(4) = i (7.i) 

i=l 

for each j. Here, as before, j is the iteration and i the translation-number index. The results will of course have to 
be averaged over the disorder ensemble, denoted by (• • •) and defined as mean values over all possible sets {a^}. Care 
must be taken to apply this averaging to the observable intensities (and not the amplitudes), since physically, one 
measures the intensities from a given realization of the disorder, and averages over the different measurements. Thus: 
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n — 1 k 



/II J_ A, 

HHdSi P(ai)\f n (q)\ 



(7.2) 



j=0 i=l 



In order to visualize the resulting random fractal, it is useful to return momentarily to the hard-wall, stepped 
surface language of Sec. IV: The support with randomized translations has steps of constant width as basic building 
blocks, but these are spaced randomly over an underlying "Cantor grid" . Due to the unequal translations, however, 
overlaps of steps may now appear, as illustrated in Figjl]. It will be shown next that in the present randomized case, 
again the intensity spectrum is sclf-amne, with the same relations between Holder exponent and fractal dimensions 
as for the non-random situation. 

The change from the deterministic case is that now the fractal operator is given by 



•^=n( i+ E^ci/ s ), 

j=i i=\ 



(7.3) 



where the random s hifts {aj} are chosen from P(a). Since one still has two identical contractions, the fractal dimension 
is unchanged [Eq. ( 3.12) )] , as required. 



A. TD Case 



Suppose a measurement is performed on a given random fractal. As for the calculations leading to the scattering 
amplitude in the deterministic case [Eq. fL7T)], the difference arises in that every translation Si is replaced by a 3 i: so 
that now: 



In 



_. n — 1 / k 



__ ? S ^< L I F{qs n ). 
j=0 \ i=l 

The resulting intensities have to be averaged over the disorder ensemble: 

2 



(7.4) 



ivrW-")(n 



\j=0 



(k- 



nd p 

^w*-) n / 



k 

1 + ^ e ls3 ^ 



i=l 



(7.5) 



For k = 1 (ID), since the shifts are chosen independently, the last expression simplifies into a product, and one obtains 
for the average intensity: 



(/«(?)> = s n I (qs n ) [] (1 + (axis' gojL)}) (ID). 
j=o 



(7.6) 



In general, no such simplification occurs, but the scaling is still tractable: 



(I n+1 (q)) = (I n (qs)) 



k+1 



1 + J2 



(7.7) 



In order to express this most accurately in the approximate general self-affine form of Eq. ( 4.12f) , the average should 
be performed at q = 0. Using the normalization condition of the distribution of translations, Eq. (7.1) one finds: 



</„+i(<7)> = s- ai (/„(?«)>, 
with a± = D — d, just as in the deterministic case [Eq. ( 6. ldj )] p8| . 



(7. 
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B. BU Case 



The scattering amplitude is now given by: 



fn(Q) 



(fe+1) 



1 n-l / k \ 
7 j=0 \ i=l / 



Averaging the intensities over the disorder ensemble: 



V ' 3=0 ^ 



The resulting scaling relation is: 

(Wj _ i W)) 

7 (g) fc + l / (5/s) 
Performing the average at q = 0, one obtains: 



n<p(a?) 



(Jn+i(g)) = g - a2 (I n (q/s)) 



W/s) 



i=l 



(7.9) 



(7.10) 



(7.11) 



(7.12) 



with «2 = £>, again as in the deterministic case [Eq. (6.15)]. 

To conclude, translational randomness alone appears to have no effect on the scaling properties of the diffraction 
spectrum. 



VIII. FURTHER PROPERTIES OF THE DIFFRACTION SPECTRUM 

The diffraction spectrum is characterized by more than just its scaling properties. Such features are discussed next. 



A. Role of Form Factor 

So far, most of the discussion has centered around the universal scaling properties of the diffraction spectrum, which 
were completely determined by the "kinematic" structure factor. However, the role of the "dynamic" form factor 
cannot be ignored in discussing the properties of the spectrum. It is in this respect that the different physical probes 
discussed in Sec.|n| differ, and that universality is broken. The form factor embodies the details of the interaction 
between probe and scatterer, and through it the potential enters the intensity spectrum. The example of He scattering 
will serve to illustrate the point. In this case, the He/surface interaction potential enters in a highly non-trivial way 
[Eq. (2.10)]. One of the striking consequences is the appearance of "rainbow" peaks in the diffraction spectrum p7[ . 
These arise essentially whenever a He atom is scattered from an inflexion point of the potential (corresponding to 
maximal force applied to the atom), typically due to an adsorbed cluster. Following is a brief discussion of the origin 
and physical significance of rainbows (see Ref. Jl9"| for a more extensive treatment). It is u sefu l to employ a stationary 
phase, approximate evaluation of the Sudden approximation scattering amplitude, Eq. ( |2.9[ ). In ID, the stationary 
phase condition is: 



q=-2r ) '(x), 

which yields x(q). The scattering amplitude is then approximated by: 

e iq-x(q) e 2iri[x(q)] 



The rainbow condition is the existence of an inflexion point in the phase shift: 

rf'{x) = 0. 



(8.1) 



(8.2) 



(8.3) 
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The poi nt X p satisfying this condition dominates the scattering by contributing a large peak. In the classical limit 
of Eq. ( |8.2| ), this shows up as a si ngularity in the intensity distribution, at momentum transfer qo satisfying the 
stationary phase condition [Eq. (^J)] together with xq. The singularity of this crude classical evaluation is smoothed 
into a finite peak in the more refined Sudden approximation calculation. 

Such Sudden approximation calculations were performed for an Ag/Pt(lll) BU Sierpinski carpet system, with a 
realistic potential, described in detail in Ref. |45|] . The results are shown in Fig.[| (fractal system) and Fig. [7] (rainbow 
analysis for a single adsorbate). It appears that, although for a single adatom the rainbow peaks are a dominant 
feature (Fig.^j), in the case of a fractal system, their role is rather negligible in determining the structure of the 
spectrum (Fig.||). The reason for this is that they are far too broad to appear as individual peaks, along with those 
due to the fractal support. The rainbows, as well as all other features of the form factor, act as very broad envelopes 
to the detailed spectral structure. The main effect of the form factor is to provide an overall intensity decrease, 
without in any way significantly altering the details of the structure factor. Since in practice one measures the full 
intensity distribution, this can have an effect on its self-affinity properties, and for a BU fractal care should be taken 
to divide by the form factor. Similar results are expected to be found in the diffraction spectra of other probes, where 
dynamical factors play an important role, but cannot lead to very peaked spectral features. 



B. "Bragg Conditions" and Band Structure 



Consider the conditions for maxima derived from the scattering amplitudes for TD and BU fractals [Eqs.(| 
For TD fractals, the condition is 



s 3 q-di 



2irti 



(6.12)] 



(8.4) 



whereas for BUs, it is 



q ■ a z 



(8.5) 



Here ij and bi are integers, and < j < n — 1. These are the "Bragg conditions" for iteratively generated fractals. 
However, since one cannot speak of a conventional unit cell with primitive lattice vectors in the fractal context, 
the pres e nt co nditions for maxima are rather different from those for periodic crystals. For ID Cantor-like sets, 
Eqs.(^4|),(^5|) reduce to: 



s 3 q ~ 



2nt 
La 

2TTb 

La 



TD 



BU 



(8.6) 



Considering first the TD case, the maxima occur for those q's which, when multiplied by s°, s , ■ ■ ■ , s" , are always 
integer multiples of j^-. For the ternary Cantor set (s = 1/3, a = 2/3), with L = 1, these g's are all the integer 
multiples of 3™7r. For 1/s equal to an arbitrary integer, these are the integer multiples of (l/s)™ -1 -^. For 1/s non- 
integer, see Ref. p7l . Cast in the usual Bragg condition language, Las 11 ^ 1 would be an effective "lattice constant" . The 
meaning of this number in the present context, is similar: it is the length of the elementary building block of the fractal 
at the n th iteration: the union of adjacent narrow black and white bars in Fig|l| (left). However, larger structures 
also repeat themselves in the fractal, with smaller frequency. These give rise to the secondary maxima in Fig.H, and 



mathematically correspond to those g's which yield integer multiples of 



271 



for only a subset of s°, s 1 



The 



incommensurability of these varying-scale, repeating structures, is what yields the multitude of peaks in the spectrum, 
as opposed to just Bragg peaks in the case of a peri o dic c rystal, and is ultimately responsible for the self-affinity of 
the spectrum. The more general conditions Eqs.( p.4[ ),(B.5), can be interpreted in a similar fashion. 

The distinction between the BU and TD cases is straightforward: the peak spacings in the former tend to zero (with 
the peak nearest to the specular found at f^s"), whereas in the latter the spacing is unbounded. The only limitation 
on the position of the furthest observable peak in the TD case is energy conservation. In both cases, however, 
the structure factors are invariant under a combination of translations and dilations (apart from the reduction in 
intensities, responsible for the self-affme properties). One is thus led to define a new basis of primitive vectors for the 
reciprocal space, from which a Brillouin zone can be constructed. As seen in Fig.|, the regions connected by these 
operations do not overlap, and can be considered as separate bands. A detailed treatment of this point is given in 
Ref. P7|, and will not be repeated here. 
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IX. SELF-AFFINE OR POWER-LAW? 



As mentioned in the Introduction, the common wisdom relating to scattering by random fractal objects (e.g. porous 
solids Jl6|,[l^]), amply confirmed experimentally, is that close to the specular the intensity satisfies a power law: 

I(q) « qT 1 (9.1) 

with 7 = £)+const. This power-law decay is clearly very different from the self-affine intensity spectrum predicted 
here for iterative fractals. Considering the unquestionable experimental evidence for the power-law, this discrepancy 
calls for clarification. The following arguments may shed some light on this issue. 

In order to derive the power-law [Eq. fl9.1| )], one typically starts with the definition of a "mass fractal dimension", 
describing the scaling of the mass N(r) enclosed in a sphere of radius r, centered at an arbitrary point in the fractal: 

N(r)^r D . (9.2) 

If the fractal is self-averaging (an assumption which is implicit in the derivation of, e.g., Refs. |l^,|l^]), then this mass 
is related to the pair distribution function g{r) by 

N(r) = {p) f g{r')d d r>, (9.3) 



with d the embedding space dimension and (p) the average density. From general scattering theory it is known that 
the structure factor is 

S(q) = 1 + (p) [ [g(r) - l]e^d d r. (9.4) 



From here one arrives at Eq. (3.1) (see Ref. |6J for details). 

It is thus seen that the crucial assumption invoked in this derivation is the scaling law Eq. (^^). It must be realized, 
however, that this expression is in many cases only true on average. This can be seen very clearly for the ternary 
Cantor set (Fig.|l|). Suppose the set has bars of unit height and one calculates its cumulative mass in the j th 
iteration, starting from the left, and in terms of the number of black bars. Then the following recursion formula may 
easily be verified: 

aKj+i) = N U) u { N U) times[last(A^)]} U {N^ + [last(A^)]}J^ 3) ' JV<°> = {1}. (9.5) 

Here liV^I is the length of the sequence N^' and last(iV^) is its last term. Indeed, the zeroth-iteration ternary 
Cantor set consists of 1 black bar, the first iteration has a cumulative mass of {1, 1, 2} black bars, the second iteration 
has mass {1, 1,2,2, 2, 2, 3, 3, 4}, etc. In Fig.(ph, is displayed on a log-log plot, together with the power-law N(r) 



[as suggested from Eq. (9.1)], i.e., a line with slope In2/ln3, the fractal dimension of the ternary Cantor set. It can 
be seen that this line serves as an accurate envelope to the actual N(r), which is in fact a "Devil's staircase", with 
a ver y ri ch (fractal) structure. This example illustrates the general situation: A simple scaling law of the form of 



Eq. (9.1) is only an average representation of the actual cumulative mass function of a fractal, whi ch m ay in fact not 
be self-averaging. Since in this work the exact properties of the fractal [i.e., equivalent to Eq. ( |9.5[ )] were used to 
calculate the scattering intensities, it should now come as no surprise that the resulting diffraction spectra themselves 



displayed the full, rich structure of the scattering fractal object. Conversely, had the power-law form of Eq. (9.1) 
been used in the present scattering calculations, the result would have been a power-law decay of the intensity. 

Why then do experiments from natural fractals yield the power-law? The preceding arguments strongly suggest 
that this is related to an averaging process which smoothes the fine-structure of the intensity distribution. A priori, 
two types of averages could be considered: (1) over the position of the center point of the cumulative-mass calculation, 
and (2) over the disorder ensemble. The first type can be ruled out immediately, however, since it is common to both 
the power-law and self-affine spectra derivations: The calculation of an intensity involves a double integral in which 
all pairs of points appear in the form J J drdr' exp[iq(r' — r)]n(r) n(r'). This automatically performs the first type 
of average. Thus, by elimination, the ensemble average appears to be responsible for smoothing out the self-affine 
properties into a simple power-law decay. Indeed, in the deterministic iterative fractals considered here, there is 
of course no ensemble to average over, in contrast to the typical experimental situation. Interestingly, the random 



fractals of Sec. VII are "not random enough", since they also display a self-affine spectrum. The type of randomness 
encountered in experiments yielding the power-law must lead, in contrast, to self-averaging between the physical 
cut-offs. There is a further difference between the randomness considered here and that encountered in experiments, 
namely that in the latter the randomness does not preserve the FD above the upper cut-off. In contrast, the type of 
randomness considered in this work preserves the FD on all scales. It is possible that this difference plays a role in 
creating the discrepancy between the experimental power-law results and the theory presented here. 
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X. CONCLUSIONS 



In summary, the diffractal- Fourier transform problem, for scattering of coherent waves from a wide class of iteratively 
constructed fractals, was solved analytically, yielding the scaling properties of the diffraction spectrum. The class 
of fractals considered here is not that which is typically observed in scattering experiments, and is characterized 
by a self-affine intensity spectrum. A simple relation was found to exist between the self-affinity exponent of this 
spectrum and the fractal dimension of the scattering fractal support. In contrast, many experiments yield intensity 
distributions characterized by a power-law decay. It is argued here that this is predominantly the result of scattering 
from self-averaging random fractals, which are more abundant in experimental realizations of fractality. The results 
apply to a large variety of scattering probes, from neutron to He scattering, the condition being the applicability 
of the Fourier transform. The differences among the probes are contained in a form factor, which, however, does 
not seem to have an important role in determining the details of the diffraction spectrum. It would be of interest 
to see whether scattering from non-self-averaging (random) fractal systems will yield a self-affine intensity spectrum 
as predicted here. Further theoretical work will concentrate on generalizing the types of randomness studied here, 
and on investigating the possible role of cut-offs in leading to the power-law decay of the intensity observed in many 
experiments. 
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FIGURE CAPTIONS 

FIG. 1. Left: Third and fourth iterations of a step on a ternary Cantor set support. Middle and right: Same, but with 
randomized translations. 

FIG. 2. Generator and first two iterations of the top-down Sierpinski carpet, supporting a harmonic potential well (contour 
lines) . By expanding each iteration so that every square is of unit size, the corresponding bottom-up fractal can be obtained. 



FIG. 3. Test of the scaling relation for TD fractals [Eq. (4.15)]: Superimposed intensities (arbitrary units), from Eq. (4.13), 



for He scattering by a hard-wall step function, built on the fifth and sixth iterations of a ternary (top) and generalized TD 
Canto r set with s = 1/5, a = 4/7 (bottom). The intensity from the fifth iteration (dotted line) is rescaled according to Eq. 



(4.15). Clearly, the rescaled intensity serves as an accurate envelope. In the n — > oo limit, therefore, subsequent iterations 
become indistinguishable and the intensity is self-affine. The insets show magnifications, in which a coarse-grained reproduction 
of the entire peak structure can be identified, illustrating the self-similarity of the spectrum. 



FI G. 4. Test of the scaling relation for BU fractals [Eq. (4.18)]: Superimposed structure factors (arbitrary units), from 
Eq. (4.11), for He scattering by a hard- wall step function , bu ilt on the same Cantor sets as in Fig.[| The intensity from the 



fifth iteration (dashed line) is rescaled according to Eq. (4.18). Again, the rescaled intensity serves as an accurate envelope, 



although the agreement worsens with increasing q. Insets as in FigM 

FIG. 5. Results of epsilon-variation analysis ^] of the intensities displayed in Fig.^J The slope of the log-log plots yields 
the self-affinity exponent as 0.633 for the n = 8 ternary Cantor set and 0.41 for the n — 6, s — 1/5, a — 4/7 set. 

FIG. 6. Top: Structure factor for He scattering from Ag adatoms on a Pt(lll) surface, with the Ag adatoms positioned 
on BU, sixth iteration Sierpinski carpet. The generator is a Pt(lll) unit-cell (L — 2.11 A). The self-similar structure can be 
noticed upon careful examination. Bottom: The complete intensity spectrum, after multiplication by the form factor (Fig.Q). 
The effect is mainly an overall intensity decrease with increasing q. The rainbows are too broad to be noticed as individual 
peaks. 
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FIG. 7. Top: Classical turning points for a single Ag atom adsorbed on a flat Pt(lll) surface, for He at normal incidence 
with k z — 6A _1 . The inflexion points are indicated (1-3), along with rays (guide for the eye only), indicating the trajectories 
of classical particles scattered from these points (note the d iffere nce in scale between the axes, causing the apparently non 



mirror-like reflection). In the hard- wall approximation [Eq. (2.12)], the inflexion points coincide with those of the phase-shift 



function, and approximately yield the positions of the rainbow peaks through the stationary phase and singularity conditions 



[Eqs.(8.1),(p,3j)]. Using this, the scattering angles are found to be (1) 21.3° and (2) 2.4° with respect to the normal to the 
surface, corresponding to q = 2.2A -1 and 0.25 h" 1 , which are approximately the rainbow positions indicated in the scattering 
intensity (bottom). 

FIG. 8. Log-log plots of the exact cumulative-mass relation [Eq. (b.5[) 1 for a seventh iteration ternary Cantor set, and a 
power law with exponent equal to this set's fractal dimension [Eq. (p.]])]. 
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